Energy Cost to Make a Hole in the Fermi Sea 
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The change in energy of an ideal Fermi gas when a local one-body potential is inserted into the 
system, or when the density is changed locally, are important quantities in condensed matter physics. 
We show that they can be rigorously bounded from below by a universal constant times the value 
given by the semiclassical approximation. 
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A problem of long-standing interest in condensed mat- 
ter physics is to give an effective estimate of the minimum 
change in kinetic energy, 5T(6p) , of an infinitely extended 
ideal Fermi gas when the density is changed locally by a 
fixed, specified amount Sp(r), i.e., the density is changed 
from a constant po to p(r) = po + Sp(r). Note that 5p(r) 
can be negative, as long as Sp(r) > — po, hence the word 
"hole" in our title. 

An equivalent problem is to calculate the minimum 
change in total energy 6£(V) of the Fermi gas when a 
local, one body potential V(r), of either sign, is added to 
the kinetic energy Hamiltonian, i.e., —V 2 — > — V 2 + V(r), 
where we use units such that h = 2m = 1. This rigorous 
equivalence, recalled in |T]) and Q below, is well known 
in density-functional theory. 

In this Letter we will give an effective answer to both 
questions by proving rigorously for dimensions D > 2 
that the well known semiclassical approximations are, 
up to an overall constant, a lower bound for the kinetic 
energy cost, as well as for the energy shift caused by 
V(r). In the semiclassical approximation one associates 
one quantum state with each box ApAr in phase space 
of volume 2tt, but this calculation is qualitatively wrong 
for D = 1 because of a singularity at the Fermi surface 
related to the Peierls instability. 

Plainly, there can not be an upper bound to 5T{Sp) be- 
cause we can always put the particles in high momentum 
states while keeping p(r) fixed. The interesting compu- 
tational question is the lower bound, i.e., the minimum 
required payment for a perturbation of po- 

The literature on S£(V) is mostly based on pertur- 
bation theoretic ideas. Such calculations are valid in 
many important cases but they do not bring the physics 
to the foreground as sharply as the semiclassical for- 
mula does. It is, therefore, important, conceptually 
and computationally, to be able to view the physically 
transparent semiclassical formula as yielding a rigorous, 
non-perturbative bound. The semiclassical formulas are 
strictly local in position-space, and hence additive over 



impurities, making them useful for density-functional 
theory [lj. No multiple-scattering calculation is needed 
here. 

The two energy shifts, ST(Sp) and S£(V), arc con- 
nected via a Legendre transform as 

ST(Sp) = sup (s£(V) - [ V(r)(po + 6p(r)) d 3 r) , (1) 

V(r) V JR 3 J 

5£(V) = inf (ST(5 P ) + / V(r)(p + 5p(v))d 3 r) . (2) 

s P( r ) V JR 3 / 

If po = then p(r) = 5p(r) > and Tip) = 5T(Sp). In 
this case we are just creating a pile of N electrons with 
density p(r) and with L 3 p(r) d 3 r — N or, equivalently, 
we are filling the negative energy states of a potential 
V. The 3D semiclassical (also known as Thomas-Fermi) 
energies are 



Tsc(p) = (3/5)(6^ 2 / 9 )s / p(r)§d 3 r, 



(3) 



£ s dV) = -(g/157r 2 ) / V(r)l d\ if Po = 0, (4) 

Jr 3 

where y± = max(0, ±y) > is the positive or negative 
part of a number y, and q is the number of available spin 
states per particle, which is 2 for unpolarized electrons. 
When po > 0, the semiclassical quantities in 3D are 

STdSp) = % c {p +5p) - % c {po) - p / <5p(r)d 3 r 

JR 3 



3 /6tt^ 
5 v q 



(p + <5p(r))* -p 3 



5 2 

-p$Sp(r) \d 3 r, 



5£ SC (V) =£ sc iV- p)-£ sc (-p) 

q 



157T 2 



(V(r)-p)l-pi 



(5) 



(6) 



where p = (67r 2 po/g) 3 is the chemical potential at den- 
sity p . 
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Main theorem. Our main Theorems in 3D are 



ST(5p) > 0.1279 8T sc (5p) 



(7) 



for the change in kinetic energy, and 



0>SS(V)-p o / V(r)d 3 r 

> 21.85 [SS sc (V)-p ( V(r)d 3 r 



(8) 



for the change in energy when an arbitrary local potential 
V(r) is inserted into the system. Similar results hold 
for all D > 2. Our bound © quantifies the validity 
of first-order perturbation theory, since po J R3 V(r) d 3 r 
is precisely the first-order term. In ([5]) note the (8p)i 
dependence for large Sp, as in ((3|), but (Sp) 2 for small S p. 

A lower bound in the pa = case was provided by the 
Lieb-Thirring inequality for all D > 1. For D = 3 

T(p) - ST(5p) > 0.6724 Tsc(p), when p = 0. (9) 

It is widely believed that 0.6724 can be replaced by 1, 
and there is continuing research in this direction. The 
inequality T(p) > KT sc (p), as in is equivalent to 



£{V) > K- 3 / 2 £ SC (V), when p Q = 0. 



(10) 



The inequality ^ was derived by first proving (|T0|) 
and then using the equivalence ((T])-([2]). Our attempt 
to follow this route in the positive density case was not 
successful. The situation changed when Rumin Q found 
a way to prove directly the kinetic energy bound ([9|) . By 
suitably modifying his method, we are now able to derive 
a lower bound for 5T{Sp), and consequently on 58 (V), 
when po > 0. 

Several papers, e.g. 0, Hj], deal with this problem from 
different points of view. In almost all cases, it has been 
approached from the side of computing the energy shift 
with a given potential V. The idea of computing the 
shift caused by a given change in density does not seem 
to have been widely considered. 

If one fixes the particle number N in a very large box 
and calculates the shift in energy caused by V, the answer 
depends on the box shape and boundary conditions 
We are able to avoid these problems and go directly to 
the thermodynamic limit by fixing the chemical potential 
p and working in infinite space. Then we can look at 
the unperturbed Hamiltonian Hq = — V 2 — p or the per- 
turbed one Hy = — V 2 — p + V(r) and fill all the negative 
energy states, i.e., all states below the Fermi level. No 
box is required in our approach; these Hamiltonians are 
defined on the whole space M. 3 . The individual energies 
are necessarily infinite but the difference is finite, as we 
will explain. 



To put this more precisely, we set Ily, respectively n , 
equal to the projections onto the negative spectrum of 
Hy and Hq. The change in energy is then 



SS(V) = Tr (H v n v - ij n 



(11) 



where Tr denotes the trace. 

To define the related kinetic energy shift, we consider 
a one-particle density matrix 7(1*, r') (suppressing spin 
indices for simplicity) and compute 



ST(Sp) =infTr H (j 

7 



no). 



(12) 



On the right side, we take the infimum over all density 
matrices 7(1*, r') whose diagonal is 7(1", r) = p + Sp(r). 
The Fermi statistics enters via the condition on 7: It is 
known Q that the necessary and sufficient condition for 
a one-body 7 to come from an iV-body fermionic den- 
sity matrix is that < 7 < 1, as an operator inequality. 
The same condition is imposed in (TT2l . Note that ST {dp) 
includes the chemical potential p in its definition (since 
H n = —V 2 — p). With this choice, ST{Sp) is always posi- 
tive, regardless of the sign of Sp(r). Formally, the reason 
for this is that IIo is the minimizer of Trf/07 among all 
one-body density matrices 7. 

Derivation of the lower bound ([7]): To simplify the 
notation, we shall assume that q = 1 from now on. 
The general case is analogous. Referring to Eq. (|12[) . 
we consider a density matrix 7(1*, r') whose density is 
Po + Sp(r) . We have to study Q = 7 — IIo, which we write 
as Q = Q+++Q—+Q+-+Q-+, where Q~ = n On , 
Q~ + = iIoQ(l — IIo), etc. In Fourier space this means 
that Q — (p,q) = Q(p 2 < p)Q(p,q)Q{q 2 < p), etc., 
with the Heaviside step function. The total change 
5p(r) of the density equals the sum of the densities of 
each of these terms, e.g., p ++ (r) = Q ++ (r,r), and so 
on. Since < 7 < 1 in the sense of operators, we have 
Q++ > and -IIo < Q~~ < 0, hence p ++ (r) > and 
— po < p (r) < 0. However, p^ (r) = p ^(r) has no 
sign a priori. 

The kinetic energy of the diagonal terms Q ±=t can be 
bounded using the method of [6j. The starting point is 
the representation 

poo 

Tr(H Q ++ ) = Tr(\H \ Q++) = / dE Tr(Q++) 

Jo 

= [ d 3 r ^ dEp+ + {v) (13) 

where Q++ = P> e Q ++ P>e, P + e + {v) - Q^ + (r,r), and 
P>e is the spectral projection of |i?o| = | — V 2 — /x| onto 
energies > E. By Schwarz's inequality and Q ++ < 1, 



< J(iP\P>eQ ++ P>eW 



(^\P<eQ++P<e\ 



< ^{^\P>eQ ++ P>e\^) + ^{^\P< e \iP) 



3 



for any ip. By taking ip to be a ^-function we obtain 
v /p++(r) < \Jp% + {t) + y/r(E), where r(.E) is the (spa- 
tial) constant density of P<e, which is easily found to 
be 

When we insert this bound on (r) into (fT3")) we obtain 
Tt(H Q ++ )> [ F(p++(r))d 3 r (14) 

JR 3 

with 



Hv) 



dE 



(15) 



The function F(y) is convex (because y (\/]y\ ~ C)\ 
is convex) and behaves like the semiclassical counterpart 
in ([5]) for small and large y. The kinetic energy of Q 
satisfies the same inequality as (fl4| . Using the convexity 
of F we obtain the bound 



Tr(ffoQ) > 2 / F 



p++(r) +p— (r) 



(16) 



For a different bound we consider the off-diagonal 



terms p -1 



P" 



Calculating in momentum space and 



using Schwarz's inequality 

(2tt)* f\p+-{r)\ 2 d 3 r = (2tt)% J p+~ (r)Q+- (r, r) d 3 r 

rf 3 g/^(p-q)0(p,q) 

^(p-q)l 2 \ 4 



d 3 p 



p 2 



< 



b 2 - /-*l 2 k 2 



(17) 



d 3 p 



q 2 >fi 



AlQ(p,q)lV-^lV-Ml 



The first square root factor on the right side can be 
rewritten as 



where 



*(fc) = | 



*(jfe) |p+-(k)| 2 d 3 fc 
d 3 p 



(18) 



(19) 



Our last task is to bound <£>(fc) from above. As a func- 
tion of k, $(fc) can be shown to be monotone decreasing. 
Thus, it attains its maximum at k = where it has the 
value $(0) = n 2 y/p. We deduce from flTTJ) that 

|p+-(r)| 2 d 3 r (20) 

\/P JR3 

< f d 3 pf AlQ(p,q)| 2 b 2 -M| 1/2 |9 2 -M| 1/2 - 



To understand the right side, we recall that Q = 7 
where < 7 < 1. Hence 

q 2 = ( 7 - n ) 2 = 7 2 - 7 n - n o7 + n 

< 7 - 7 n - n o7 + n = Q++ - c 

and, therefore, 

Tr H a Q > Tr \H \Q 2 > 2 x (right side of flU}). 
We deduce that 

16tt 



no 



TvHoQ > 



^P 



\p + -(r)\U 3 r 



(21) 
1). 



So far we have found two lower bounds, (|16l) and 
for the shift in kinetic energy ST(Sp) = inf Tr HoQ, which 
we average with coefficients t and 1 — t. Thus, ST(5p) is 
bounded below by 



2tF 



p++(r)+p--(r) 
2 



ig (l-i)|p+-(r)| 2 



where _F is in (1151) . We must now give a lower bound 
to the right side of (|22|) in terms of the total change 
in density 5p(r) = p ++ (r) + p (r) + 2p~ + (r). These 
three quantities are not known separately but they do 
satisfy the constraints that p ++ (r) + p~~(r) > — po and 
<5p(r) > ™ Po- We then look for a < t < 1 such that 

16tt 



2tF 

\2, 
3 /6^ 2 

> ft- 



/x\ Ion . , 9 



(po + x + 2y) 3 - pi - -pi (x + 2y) 

holds for all x + 2y > —po and all x > —po, and with 
k as large as possible. Solving this problem numerically 
leads to t = 0.7267 and k = 0.12797. This completes the 
derivation of our first main result ([7]). 

Using ([2]) we obtain the bound on the shift in energy 
S£(V) > k S£ sc (ft -1 ^). Finally we can use the fact 
that S£ SC (V) — po J* R3 V(r)d 3 r is a monotone decreasing 
function of p. This implies (|2"3"|) . which implies ©: 

-El [ y (r)d 3 r 

k- 5/2 U£ SC (V) - po V(t) d 3 r^j . (23) 



S£ sc I -V 
ft 

> 



Extension to 2D: Our method can be generalized to 
2D (indeed to any dimension except ID). The result is 

6T 2D {Sp)> 0.04493 6T 2D (Sp) 



(24) 



for the change in kinetic energy, and 

> 8£ 2U (V) -po / V(r)d 2 r 

Jr 2 

> 22.25 (5£ 2 ?(V)-p f U(r)d 2 r' 

V JR 2 J 



(25) 
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for the change in energy when a potential V(r) is in- 
serted into the system. The 2D semiclassical tunc- 
tions are T s 2 c D (p) = 5T™(p) = (2^/q) J R2 p(r) 2 d 2 r and 
5£™(V) = -q/(8n) J R2 ((V(r) - p) 2 _ - p 2 ) d 2 r with 
|U = Airpo/q. 

Peierls Instability in ID: In a ID free Fermi g as, a 
bound like © cannot hold, in general. When a potential 
V(r) is inserted into the system at positive density po, 
the corresponding variation of the semiclassical energy is 



6£™(V) = 



2q 
3tt 



((V(r)-pf 2 -p^ 2 )dr 



with p = (irpo/q) 2 . For small V(r), this gives 



6£™(V)-po V(r)dm- 



4ir 2 p 



V{rf dr. 



On the other hand, second-order perturbation theory 
flo| predicts that the true shift S£ 1D (V), for small V, is 



8£ m (V)-p l V(r)dr 



-2- dk 



df: 



q_ 

47T 



The logarithm diverges at \k\ — 2*//I; hence the second- 
order term can be made arbitrarily large while keeping 
J V(r) 2 dr fixed. Thus, there cannot be any lower bound 
involving S£^(V). This divergence in ID is well known, 
and is related to the Peierls instability [lo| . In higher 
dimensions, the second-order approximation is bounded 
(this follows from our bound ©), but it is known to 
have an infinite derivative at |fc| = 2y r p, a fact that is 
sometimes called the Migdal-Kohn anomaly 
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Extension to positive temperature: The change in 
free energy at temperature T = (fce/3) -1 and chemical 
potential p is 



5T(V) = -jS _1 Tr [ln(l 



ln(l 



)] 



with Hy and Hq as before. Using the fact that ln(l 



-PE\ _ 



ST(V) = /3 



(A - E) dX we find 



( 1 + e M)2 



5£» +x {V) dX 



with 5£^ +x (V) = - Tr [{H v - A)_ - (H - A)_]. This 
formula expresses the positive temperature energy shift 
as a mixture of zero temperature energy shifts with dif- 
ferent chemical potentials. In fact, S£ fJ,+x (V) is nothing 
but the energy shift estimated before, with chemical po- 



tential p + X instead of p. Thus ([8]) leads to 
> 5F{V) - pr I V(r)d 3 r 

JS.3 

>_(21.85)-|^ [ dX [ d 3 r(l + e^)- 2 e^ 

x {(V(r) -p-X) i -(p + X)l + l(p + A)Ju(r)} 
with density 



Pt = 



(27T) 



(1 + e ^(P 2 -P))-i d ; 



Similar results hold for all D > 2. This is a bound on 
the change in free energy after insertion of a potential V; 
the corresponding version in terms of the density change 
bp can be obtained via a Legendre transform, as in ([1]). 

Extension to periodic background potentials: Our 

method also works here. The result depends on knowing 
two things: the density of states close to the Fermi level, 
p, and the non- homogeneous background density po(r) 
for this p. With these quantities in hand, the calculation 
follows along the same lines as the one given here. There 
are various possible energy-band scenarios and, for lack 
of space, we defer the details to a forthcoming paper. 

Conclusion: We show rigorously that the energy shift of 
a Fermi gas, caused either by a local density perturbation 
or by a local potential, is described, qualitatively, by a 
semiclassical calculation. 
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